The impact of non-Gaussianity when searching for Primordial Black Holes with LISA

This paper demonstrates that the presence of astrophysical foregrounds and strong primordial non-Gaussianity significantly weakens the link between scalar-induced gravitational waves and primordial black hole abundance, thereby limiting LISA's ability to definitively confirm or rule out asteroid-mass primordial black holes as dark matter.

The Gist

The Big Picture: Listening to the Baby Universe

Imagine the Universe as a giant, quiet room. For a long time, we've been trying to hear the faint whispers of the very first moments after the Big Bang. One of the best ways to do this is by listening for Gravitational Waves (GWs)—ripples in the fabric of space-time, like sound waves traveling through water.

The paper focuses on a specific space telescope called LISA (Laser Interferometer Space Antenna), which is like a giant, floating ear in space designed to listen to a specific "pitch" of these ripples (the millihertz range).

The scientists are asking a big question: Can LISA tell us if the Universe is full of tiny, ancient black holes (called Primordial Black Holes or PBHs) that make up Dark Matter?

The Two Clues: The Echo and the Shadow

The paper argues that there are two clues left behind by the same event in the early Universe:

1. The Echo (SIGW): When the early Universe was turbulent, it created "Scalar-Induced Gravitational Waves" (SIGWs). Think of this as the echo of a loud shout. If the early Universe was chaotic enough to make black holes, it must have also made this echo.
2. The Shadow (PBHs): These are the Primordial Black Holes themselves. If the early Universe had big enough "bumps" in density, they would have collapsed into black holes.

The Standard Logic:
In the past, scientists thought the relationship was simple: If you hear the Echo (SIGW), you know the Shadow (PBHs) exists. If you don't hear the Echo, the Shadows don't exist.
If LISA listens and hears nothing, the old theory says: "Okay, the asteroid-mass black holes that could be Dark Matter definitely aren't there."

The Twist: The "Weather" of the Universe (Non-Gaussianity)

Here is where the paper introduces a game-changer. The scientists realized that the early Universe might not have been "smooth" or "average" (Gaussian). It might have had weird, extreme weather patterns. In physics, this is called Non-Gaussianity (or $f_{NL}$).

The Analogy: The Lottery Ticket
Imagine the early Universe is a lottery machine that prints numbers (density bumps).

• Gaussian (Normal): The machine prints numbers in a perfect bell curve. Most numbers are average; extreme numbers are incredibly rare.
• Non-Gaussian (Weird): The machine is rigged. It might print mostly average numbers, but occasionally, it prints a "Jackpot" number that is massive.

The paper shows that Non-Gaussianity changes the rules of the lottery:

• The Echo (SIGW): The echo depends on the average size of the numbers. It doesn't care much if there are a few Jackpots. It's like hearing the general hum of a crowd; a few people shouting doesn't change the overall noise level much.
• The Shadow (PBHs): The black holes only form from the Jackpots (the extreme outliers). If the machine is rigged to produce more Jackpots (positive Non-Gaussianity), you get way more black holes, even if the average noise (the Echo) stays the same.

The Main Discovery: The "Loophole"

The authors found that Non-Gaussianity creates a massive loophole.

1. The Old View: If LISA doesn't hear the Echo, we can rule out the black holes.
2. The New View: If the early Universe was "rigged" (had high Non-Gaussianity), the Echo might be too quiet for LISA to hear, BUT the black holes could still be everywhere because the "Jackpots" were so extreme.

In simple terms: You can have a room full of black holes (Shadows) without making a loud enough echo (SIGW) for LISA to hear, if the distribution of the black holes is weird enough.

The Problem: The "Blind Spot"

The paper also highlights a frustrating problem. Even if LISA does hear the Echo and measures it perfectly, we still can't be sure how many black holes there are.

The Analogy: The Volume Knob
Imagine you are trying to guess how many people are in a stadium based on the noise level.

• If you know the "volume knob" (Non-Gaussianity) is set to 0, you can guess the crowd size accurately.
• But if you don't know where the volume knob is set, the same noise level could mean a crowd of 10 people or a crowd of 10 billion people.

The paper calculates that because we don't know the exact value of Non-Gaussianity ($f_{NL}$), our estimate of the number of black holes could be off by 30 orders of magnitude. That's the difference between a few grains of sand and the entire Earth!

What Does This Mean for LISA?

1. LISA is still a great detective: It will be able to measure the "Echo" (the gravitational waves) very precisely.
2. But it can't solve the mystery alone: Without knowing the "weather" of the early Universe (the Non-Gaussianity), LISA cannot definitively say "Yes, black holes are Dark Matter" or "No, they aren't."
3. The "Asteroid Window": There is a specific size range of black holes (asteroid-sized) that LISA is perfect for finding. The paper shows that if we don't find the Echo, we can't rule out these black holes anymore because the "Non-Gaussianity loophole" keeps them hidden.

The Conclusion

The paper is a warning and a guide. It tells us that we cannot ignore the "weirdness" (Non-Gaussianity) of the early Universe.

If we want to use LISA to prove or disprove that tiny black holes are Dark Matter, we need to understand the "shape" of the early Universe's fluctuations much better. Otherwise, LISA might listen to the silence, and we might wrongly conclude the black holes aren't there, when in reality, they are just hiding behind a curtain of non-Gaussian statistics.

Summary in one sentence:
LISA is a powerful microphone that can hear the echoes of the early Universe, but because the early Universe might have had "weird weather" (Non-Gaussianity), hearing silence doesn't necessarily mean the black holes aren't there—they might just be hiding in the noise.

Technical Summary

1. Problem Statement

The paper addresses the potential of the Laser Interferometer Space Antenna (LISA) to constrain the abundance of Asteroid-Mass Primordial Black Holes (PBHs) as a dark matter candidate.

• The Link: In standard cosmological scenarios, large primordial curvature perturbations ($\zeta$) on small scales lead to two observable consequences:
  1. Primordial Black Hole (PBH) formation via gravitational collapse of overdensities.
  2. Scalar-Induced Gravitational Waves (SIGWs) generated at second order in perturbation theory.
• The Standard Assumption: Previous studies often assume primordial perturbations are Gaussian. Under this assumption, the relationship between the amplitude of the scalar power spectrum ($A$) and the resulting PBH abundance ($f_{\text{PBH}}$) is well-defined. Consequently, a non-detection of SIGWs by LISA was thought to definitively exclude the asteroid-mass PBH dark matter window.
• The Gap: This conclusion relies heavily on the assumption of Gaussianity. The paper investigates how primordial non-Gaussianity (NG), parametrized by $f_{\text{NL}}$, and astrophysical foregrounds alter the link between SIGWs and PBHs, potentially invalidating the "non-detection = exclusion" conclusion.

2. Methodology

The authors employ a multi-faceted approach combining analytical calculations, numerical simulations, and Bayesian inference:

• Theoretical Framework:

  • Non-Gaussianity Model: They adopt the local type of non-Gaussianity, expressing the curvature perturbation $\zeta$ as a function of an auxiliary Gaussian field $\zeta_G$: $\zeta = \zeta_G + \frac{3}{5}f_{\text{NL}}(\zeta_G^2 - \langle \zeta_G^2 \rangle)$.
  • Power Spectrum: They assume a log-normal primordial power spectrum characterized by amplitude $A$, peak wavenumber $k_*$, and width $\Delta$.
  • SIGW Calculation: They calculate the SIGW energy density ($\Omega_{\text{GW}}$) by decomposing the 4-point correlation function of $\zeta$ into disconnected (Gaussian) and connected (NG) components. The NG contribution scales as $f_{\text{NL}}^2$.
  • PBH Calculation: They use Threshold Statistics (TS) to estimate the PBH abundance ($f_{\text{PBH}}$), integrating the probability distribution of density contrasts above a critical threshold. They also briefly compare with Peaks Theory (PT).
  • Perturbativity Check: They rigorously define the validity of the perturbative expansion, establishing a condition $A f_{\text{NL}}^2 \lesssim 0.1$ to ensure the local expansion remains valid.

• LISA Sensitivity Analysis:

  • They calculate the Signal-to-Noise Ratio (SNR) for LISA over a 4-year observation, incorporating both instrumental noise and astrophysical foregrounds (Galactic white dwarf binaries and extra-galactic compact binaries).
  • They derive projected constraints on the scalar power spectrum amplitude $A$ for various $f_{\text{NL}}$ values and spectral widths.

• Bayesian Inference:

  • Using the SIGWAY and SGWBinner pipelines, they generate mock data for injected SIGW signals.
  • They perform Markov Chain Monte Carlo (MCMC) scans to reconstruct the posterior distributions of the power spectrum parameters ($A, \Delta, k_*, f_{\text{NL}}$) and subsequently infer the derived PBH quantities ($f_{\text{PBH}}, \langle M_{\text{PBH}} \rangle$).

3. Key Contributions

1. Decoupling of SIGW and PBH Constraints: The paper demonstrates that while SIGW detection limits depend primarily on $|f_{\text{NL}}|$ (quadratic dependence), PBH formation is exponentially sensitive to the sign and magnitude of $f_{\text{NL}}$.

   • Positive $f_{\text{NL}}$: Enhances the tail of the probability distribution, drastically increasing PBH abundance for a fixed $A$. This allows PBHs to constitute dark matter even if the SIGW signal is below LISA's detection threshold.
   • Negative $f_{\text{NL}}$: Suppresses the tail, reducing PBH abundance.

2. Re-evaluation of the "Exclusion Window":

   • In the Gaussian limit ($f_{\text{NL}}=0$), a non-detection of SIGWs by LISA would exclude PBHs in the asteroid-mass window ($10^{-17} M_\odot \lesssim M_{\text{PBH}} \lesssim 10^{-5} M_\odot$) as dark matter.
   • With large positive non-Gaussianity ($f_{\text{NL}} \gtrsim 50$), the required amplitude $A$ to produce $f_{\text{PBH}}=1$ drops significantly. Consequently, LISA could fail to detect SIGWs even if PBHs make up 100% of dark matter, effectively reopening the asteroid-mass window.

3. Quantification of Uncertainty:

   • The authors show that even if LISA detects a SIGW signal with high precision (pinpointing $A, \Delta, k_*$ to sub-percent accuracy), the inference of $f_{\text{PBH}}$ remains highly uncertain.
   • Because $f_{\text{PBH}}$ depends exponentially on $f_{\text{NL}}$, and LISA's ability to constrain $f_{\text{NL}}$ is limited (especially for low SNR or specific spectral shapes), the uncertainty in $f_{\text{PBH}}$ can span 30 orders of magnitude.

4. Perturbativity Limits: The study provides a detailed map of where the perturbative treatment of non-Gaussianity breaks down ($A f_{\text{NL}}^2 \gtrsim 0.1$), highlighting regions where theoretical predictions for SIGWs and PBHs become unreliable.

4. Key Results

• Sensitivity Curves: LISA can constrain the scalar power spectrum amplitude $A$ down to $\sim 10^{-4}$ (for narrow spectra) in the Gaussian case. However, the presence of large $f_{\text{NL}}$ shifts the sensitivity curves only slightly for the SIGW signal itself but drastically changes the $f_{\text{PBH}}=1$ exclusion lines.
• Foreground Impact: Astrophysical foregrounds moderately weaken LISA's sensitivity to SIGWs but have a negligible effect compared to the impact of non-Gaussianity on the PBH abundance inference.
• MCMC Reconstruction:
  • Mean Mass ($\langle M_{\text{PBH}} \rangle$): Can be reconstructed reliably (narrow posterior) as it depends mainly on the spectral peak and width.
  • Abundance ($f_{\text{PBH}}$): The posterior is extremely broad. For an injected signal with $f_{\text{NL}}=0$, the 95% confidence interval for $f_{\text{PBH}}$ can range from negligible to $f_{\text{PBH}} \approx 1$ simply due to the uncertainty in $f_{\text{NL}}$.
• The "Incredulity Limit": The paper identifies a region where the expected number of PBHs in the Hubble volume is less than one, rendering them observationally irrelevant. Non-Gaussianity significantly shifts this boundary.

5. Significance and Conclusions

• Paradigm Shift: The paper fundamentally alters the interpretation of future LISA data regarding PBH dark matter. It argues that LISA alone cannot definitively rule out the asteroid-mass PBH dark matter scenario without independent constraints on primordial non-Gaussianity.
• Critical Role of $f_{\text{NL}}$: The uncertainty in $f_{\text{NL}}$ is the dominant systematic error in translating SIGW limits into PBH constraints. A non-detection of SIGWs does not necessarily mean "no PBHs"; it could mean "PBHs exist but are hidden by non-Gaussianity."
• Future Directions: The authors conclude that while LISA will significantly narrow the viable parameter space, a robust confirmation or exclusion of PBH dark matter requires:
  1. Better constraints on $f_{\text{NL}}$ from other probes (e.g., CMB, large-scale structure, or future GW detectors).
  2. A deeper understanding of the perturbative limits of non-Gaussian models.
  3. Joint analysis of SIGW signals with other cosmological observables.

In summary, the paper serves as a crucial cautionary note for the interpretation of future gravitational wave data, emphasizing that non-Gaussianity acts as a "degeneracy breaker" that can hide PBH dark matter from SIGW detection limits.