Weak lensing of bright standard sirens: prospects for $σ_8$

This paper demonstrates that incorporating weak lensing into bright standard siren analyses enables the measurement of the matter perturbation parameter $\sigma_8$, with future observatories like the Einstein Telescope and LISA potentially achieving 10% and 30% accuracy, respectively, given sufficient populations of neutron star and massive black hole binaries with electromagnetic counterparts.

The Gist

Here is an explanation of the paper using simple language and creative analogies.

The Big Idea: Listening to the Universe's "Echoes"

Imagine the universe is a giant, dark concert hall. For decades, astronomers have been trying to figure out how big the hall is and how fast it's expanding by looking at the lights on the stage (stars and supernovas). This is like using a "cosmic ruler" made of light.

But recently, we've started listening to the music instead: Gravitational Waves. These are ripples in space-time caused by massive events, like two black holes or neutron stars crashing into each other. Because we know how loud these crashes should be (based on physics), we can tell how far away they are just by how quiet they sound when they reach us. These are called "Standard Sirens" (like "Standard Candles," but for sound).

The Problem: The "Foggy Window" Effect

Here's the catch: As these sound waves travel across the universe to reach us, they pass through a cosmic landscape filled with invisible mountains and valleys made of dark matter and galaxies.

Think of the universe as a foggy window. Sometimes the glass is thick, and sometimes it's thin.

• If the sound passes through a "thick" patch (a cluster of galaxies), it gets magnified (louder).
• If it passes through a "thin" patch (empty space), it gets demagnified (quieter).

This is called Weak Lensing. It makes the distance measurements slightly "fuzzy." Usually, astronomers try to fix this fuzziness to get a clearer picture of the universe's expansion.

The Twist: Using the Fuzziness as a Tool

This paper says: "Wait a minute! Don't just try to fix the fuzziness. Let's study it!"

The author, Ville Vaskonen, suggests that the pattern of this fuzziness actually tells us something new. It's like looking at the ripples on a pond. If you know how the water is moving, you can figure out how many fish are swimming underneath, even if you can't see the fish directly.

By analyzing how much the "Standard Sirens" get magnified or demagnified, we can measure $\sigma_8$ (sigma-8).

• What is $\sigma_8$? Think of it as the "Clumpiness Meter" of the universe. It tells us how much matter is bunched together in big lumps (like galaxy clusters) versus how much is spread out evenly.
• Why does it matter? Standard sirens usually tell us how fast the universe is expanding. But they don't tell us much about how the stuff inside the universe is arranged. This new method uses the "fuzziness" to measure the arrangement of the cosmic web.

The Plan: Two Different Microphones

The paper looks at two future "microphones" (detectors) that will listen to these cosmic crashes:

1. ET (Einstein Telescope): A super-sensitive detector on Earth.

   • The Target: 300 collisions of Neutron Stars (the dense cores of dead stars).
   • The Result: If we catch 300 of these, we can measure the "Clumpiness" ($\sigma_8$) with 10% accuracy. That's like guessing the weight of a person within 10 pounds.

2. LISA: A space-based detector (a giant triangle of satellites).

   • The Target: 12 collisions of Massive Black Holes (the giants at the centers of galaxies).
   • The Result: Even with just 12 of these huge events, LISA can measure the "Clumpiness" with 30% accuracy. It's a smaller sample, but the events are so massive and loud that the signal is very clear.

The Analogy: The Rainstorm

Imagine you are standing in a field during a rainstorm, trying to guess how many raindrops are falling.

• Old Method: You count the drops hitting a bucket. This tells you the total rain (the expansion of the universe).
• New Method: You notice that some drops are hitting your hat harder than others because of the wind (the cosmic structures). By studying how the wind is pushing the rain around, you can figure out how many trees and bushes are in the field (the clumpiness of matter), even if you can't see them.

Why This is a Big Deal

1. New Superpower: It turns gravitational waves from just a "ruler" into a "map maker." We can now map the invisible scaffolding of the universe.
2. Cross-Checking: We currently measure the "Clumpiness" using light (telescopes looking at galaxies). This method uses sound (gravitational waves). If the two methods agree, we are on the right track. If they disagree, it might mean our understanding of gravity or dark matter is wrong!
3. Feasibility: The paper proves that we don't need millions of events to do this. Just a few hundred (or even a dozen for the biggest black holes) are enough to get a good measurement.

The Bottom Line

This paper is a proposal to stop treating the "noise" in our gravitational wave signals as a mistake to be fixed. Instead, we should treat that noise as a secret message about how the universe is structured. By listening carefully to how the universe "warps" the sound of crashing stars, we can finally get a better handle on how clumpy the cosmos really is.

Technical Summary

Here is a detailed technical summary of the paper "Weak lensing of bright standard sirens: prospects for $\sigma_8$" by Ville Vaskonen.

1. Problem Statement

Gravitational wave (GW) events with electromagnetic (EM) counterparts, known as bright standard sirens, provide direct measurements of luminosity distances ($D_L$) independent of the cosmic distance ladder. However, as GW detectors (like the Einstein Telescope, ET, and LISA) become more sensitive, the propagation of GWs through intervening cosmic structures introduces weak gravitational lensing.

This lensing causes magnification or demagnification of the GW signal, introducing a scatter in the observed luminosity distances. Traditionally, this scatter is treated as a source of noise or systematic error that must be marginalized over to measure the mean Hubble diagram (parameters $H_0$ and $\Omega_M$).

The Core Problem: The paper argues that this scatter is not merely noise but encodes valuable information about the distribution of matter in the universe. Specifically, the variance and non-Gaussianity of the magnification distribution depend on the standard deviation of matter perturbations, $\sigma_8$. Previous studies often relied on linear convergence power spectra, which are valid only on large scales and require massive sample sizes. This work seeks to exploit the full non-Gaussian magnification probability distribution function (PDF), including non-linear structures, to constrain $\sigma_8$ with smaller, realistic sample sizes.

2. Methodology

A. Modeling the Lensing Magnification Distribution ($dP/d\mu$)

The authors developed a stochastic approach (implemented in C++) to compute the magnification PDF, $dP/d\mu$, rather than relying on ray-tracing through N-body simulations (which are computationally expensive).

• Line-of-Sight Modeling: The matter distribution is modeled as an ensemble of discrete lenses (halos and filaments) drawn from a specified mass function.
• Lens Components:
  • Halos: Modeled using pseudo-elliptical Navarro-Frenk-White (NFW) profiles. Parameters include mass ($M$), redshift ($z$), projected distance ($r$), and ellipticity orientation.
  • Filaments: Modeled as cylindrical structures with a specific radial density profile ($\rho \propto (1 + r^2/r_s^2)^{-1}$).
• Mass Function: The number density of lenses is derived using the excursion-set formalism with mass-dependent collapse barriers. The model includes a linear bias factor to account for halo clustering.
• Calculation: For a source at redshift $z_s$, the authors generate realizations of lenses along the line of sight. They sum the convergence ($\kappa$) and shear ($\gamma$) contributions from individual lenses to calculate the total magnification $\mu = [(1-\kappa)^2 - \gamma^2]^{-1}$.
• Validation: The resulting variance of luminosity distances ($\sigma_{D_L}/D_L$) as a function of redshift was compared against results from high-resolution N-body simulations (Takahashi et al. 2011), showing reasonable agreement, particularly when filaments and ellipticity are included.

B. Hubble Diagram Reconstruction

The authors constructed benchmark catalogues to simulate observations by future detectors:

• ET (Einstein Telescope): 300 Binary Neutron Star (BNS) mergers with EM counterparts ($z < 2$).
• LISA: 12 Binary Massive Black Hole (BMBH) mergers with EM counterparts ($z < 10$).
• Redshifts: Assumed to be known precisely from EM counterparts.
• Distance Uncertainties: Intrinsic measurement errors were set to 3% for BNS and 0.3% for BMBH.

Inference Framework:
A Markov Chain Monte Carlo (MCMC) analysis was performed to infer cosmological parameters $\theta = \{h, \Omega_M, \sigma_8\}$.

• Likelihood: The likelihood function incorporates the full non-Gaussian magnification distribution $p_\mu(\mu|z)$ rather than a Gaussian approximation.
• Model: The probability of observing a distance $D_L$ at redshift $z$ is derived by convolving the true distance-redshift relation with the lensing magnification PDF.
• Priors: Uniform priors were used for $h$, $\Omega_M$, and $\sigma_8$.

3. Key Contributions

1. First Implementation of Full Magnification Distribution: This is the first study to implement the full non-Gaussian magnification distribution (including non-linear structures like filaments and elliptical halos) directly into the Hubble diagram reconstruction for standard sirens.
2. $\sigma_8$ Constraints from Small Samples: The paper demonstrates that $\sigma_8$ can be constrained using the scatter of the Hubble diagram, not just the mean.
3. Complementarity: It establishes that weak lensing transforms bright sirens from purely geometric probes into tracers of cosmic structure, offering a probe complementary to CMB, galaxy clustering, and traditional weak lensing surveys.

4. Key Results

Measurement Precision

• Einstein Telescope (ET): With a sample of 300 BNS systems, ET can measure $\sigma_8$ with an accuracy of 10%.
• LISA: With a sample of only 12 BMBH systems, LISA can measure $\sigma_8$ with an accuracy of 30%.
• Combined: A combination of ET and LISA data could achieve an accuracy of 8%.

Parameter Degeneracies

• $\sigma_8$ vs. $h, \Omega_M$: The analysis reveals that $\sigma_8$ has much weaker correlations with the Hubble constant ($h$) and matter density ($\Omega_M$) compared to the correlations between $h$ and $\Omega_M$. This allows $\sigma_8$ to be constrained independently of the mean distance-redshift relation.
• Bias Check: The inference was found to be unbiased; the posterior distributions centered correctly on the fiducial values used to generate the mock catalogues.

Impact of Modeling

The study highlights that including filaments and elliptical halos significantly affects the predicted weak lensing scatter. Simplified models (e.g., spherical halos only) underestimate the scatter, potentially leading to biased parameter estimates.

5. Significance and Future Outlook

• New Cosmological Probe: This work opens a new channel for studying cosmological perturbations using GWs. It allows for the measurement of matter fluctuations at low-to-intermediate redshifts, a regime where CMB constraints are indirect.
• Robustness: While the current modeling of filaments and halos is relatively simplistic, the results demonstrate the potential of the method. Future refinements in modeling non-linear structures will be crucial for reaching the precision of CMB constraints.
• Beyond $\Lambda$CDM: The approach is highly sensitive to the nature of dark matter. For instance, "warm" or "fuzzy" dark matter models suppress small-scale structures, which would reduce the lensing scatter, whereas compact dark matter objects would enhance the high-magnification tail. Thus, bright sirens could provide novel constraints on dark matter properties.
• Methodological Advancement: The paper suggests that machine learning could accelerate the inference pipeline, allowing for systematic exploration of how uncertainties scale with sample size and redshift distribution.

In conclusion, the paper establishes that weak lensing scatter in bright standard sirens is a powerful, previously underutilized tool for measuring $\sigma_8$, offering competitive constraints with relatively small sample sizes from next-generation GW observatories.