Binary Black Holes population synthesis based on the current LVK observations

This paper analyzes LIGO-Virgo-KAGRA observations from 2015 to early 2024 to demonstrate that the binary black hole population is best explained by a combination of first-generation stellar remnants, hierarchical merger products, and potentially primordial black holes, rather than solely by the collapse of massive stars.

The Gist

Imagine the universe as a giant, cosmic dance floor. For the last decade, scientists have been listening to the "music" of this dance floor using giant ears called LIGO, Virgo, and KAGRA. These ears detect ripples in space-time called gravitational waves, which are created when two massive objects, like black holes, crash into each other.

This paper is like a detective story. The authors, Mehdi, Ilias, and Muhsin, are trying to figure out who is dancing on this floor and where they came from. They have a list of over 150 "dance partners" (binary black holes) that have been caught in the act.

Here is the breakdown of their investigation, explained simply:

The Three Suspects (The Populations)

The scientists propose that the black holes they see come from three different "families" or origins. Think of them as three different types of dancers:

1. The "First-Gen" Dancers (Stellar Black Holes):

   • Who they are: These are the "normal" black holes. They are born when a massive star (like a giant sun) runs out of fuel, collapses, and explodes.
   • The Analogy: Imagine a bakery. Every time a baker makes a loaf of bread, it's slightly different, but they generally follow a standard recipe. Most loaves are small, some are medium, and very few are huge. The scientists expected these black holes to follow a similar "recipe" (a power-law distribution), where smaller ones are common and huge ones are rare.
   • The Problem: The data shows some black holes are way too heavy to be made by this standard recipe. It's like finding a loaf of bread the size of a house in a normal bakery.

2. The "Re-Merger" Dancers (Hierarchical Mergers):

   • Who they are: These are black holes that have already danced once before. Two "First-Gen" black holes crash, merge into one bigger black hole, and then that new giant merges again.
   • The Analogy: Imagine a game of "Musical Chairs" but with black holes. Two black holes merge to make a bigger one. That bigger one gets pushed into another collision. It's a "black hole family tree" where the parents are already the result of a marriage.
   • Why it matters: This explains the "giant loaves of bread." If you keep merging black holes, you get super-massive ones that the standard bakery recipe can't explain.

3. The "Ancient" Dancers (Primordial Black Holes - PBHs):

   • Who they are: These are the "ghosts." They didn't come from dead stars at all. They were formed in the very first split-second of the Big Bang, like dust motes frozen in time.
   • The Analogy: Imagine the universe is a giant ocean. The "First-Gen" black holes are like icebergs that formed from the water freezing later. The "Ancient" black holes are like rocks that were thrown into the ocean before the water even existed. They are made of "dark matter" (the invisible stuff holding galaxies together).
   • The Clue: If these exist, they would have a very specific weight and would be dancing in a very specific way that differs from the star-born ones.

The Investigation (The Method)

The authors took the list of 150+ black hole collisions and ran them through a massive computer simulation. They asked: "Which mix of these three suspects best explains the data we actually heard?"

They used a statistical tool (like a very sophisticated scale) to weigh different theories:

• Theory A: It's just the "First-Gen" bakery.
• Theory B: It's the bakery plus the "Re-Merger" family tree.
• Theory C: It's the bakery, the family tree, and the "Ancient" ghosts.

The Verdict (The Results)

The data spoke loud and clear. Here is what they found:

1. The Bakery Alone Isn't Enough: The simple "First-Gen" theory (just stars dying) failed to explain the heavy black holes. The math just didn't add up.
2. The "Re-Merger" Family is Real: The data strongly suggests that a significant chunk of these black holes (maybe 10% to 50% of the total) are "second-generation" or even "third-generation" black holes. They are the result of previous crashes. This explains why we see black holes weighing 30 to 45 times the mass of our Sun—too heavy for a single star to make, but perfect for a black hole that has already merged once.
3. The "Ancient" Ghosts are Likely: The best fit for the data came when they added the Primordial Black Holes (PBHs) to the mix.
   • The authors found that a small percentage (a few percent) of the black holes we see might be these ancient, dark-matter-born objects.
   • If this is true, it means that 0.3% to 2% of all the "dark matter" in the universe is made of these stellar-mass black holes.

The Big Picture

Think of the universe's black hole population like a fruit salad.

• For a long time, we thought the salad was just apples (stars dying).
• But when we tasted the salad, we realized there were also oranges (black holes that merged with other black holes).
• Now, this paper suggests there might even be a few exotic berries (primordial black holes) mixed in that we didn't know were there.

Why does this matter?
If the "exotic berries" (Primordial Black Holes) are real, it changes our understanding of the universe's history. It would mean that dark matter isn't just some invisible fog, but is actually made of tiny, invisible black holes formed at the birth of time. It also tells us that black holes are "social," constantly crashing and merging to build even bigger monsters.

In short: The universe is a chaotic dance floor where black holes are born from dying stars, but they also keep crashing into each other to make giants, and there's a good chance some of the dancers are ancient ghosts from the very beginning of time.

Technical Summary

Here is a detailed technical summary of the paper "Binary Black Holes population synthesis based on the current LVK observations" by El Bouhaddouti, Cholis, and Aljaf.

1. Problem Statement

The LIGO-Virgo-KAGRA (LVK) collaboration has detected over 150 binary black hole (BBH) merger events. A central challenge in gravitational-wave astrophysics is determining the formation channels of these binaries. Specifically, the authors aim to distinguish between three potential populations contributing to the observed merger rate:

1. First-Generation (1G) Black Holes: Formed from the core collapse of massive stars, expected to follow a power-law mass distribution with a potential cutoff due to the pair-instability supernova (PISN) mechanism.
2. Hierarchical Mergers (HM): Black holes formed from the merger of previous black hole binaries (2nd, 3rd, or higher generations), typically occurring in dense stellar environments like globular clusters. These are expected to produce more massive black holes.
3. Primordial Black Holes (PBHs): Black holes formed in the early universe from density perturbations, potentially contributing to dark matter and merging within dark matter halos or the intergalactic medium.

The paper investigates whether the current LVK data (O1–O4a) can be explained solely by 1G stars or if statistical evidence requires the inclusion of HM and/or PBH populations.

2. Methodology

Data Set

• Source: The authors utilize the GWTC-4.0 catalog (and related O1–O4a data), selecting 159 BBH events with a signal-to-noise ratio (SNR) > 8 and secondary mass $m_2 > 4 M_\odot$.
• Data Processing: They construct histograms for primary mass ($m_1$), secondary mass ($m_2$), and redshift ($z$) by stacking the posterior probability density functions (PDFs) provided by the LVK collaboration.

Population Models

The authors model the merger rates and mass distributions for three distinct populations:

1. First-Generation (1G) Models:
   • PL (Power-Law): $dN/dm_1 \propto m_1^{-\alpha}$.
   • PLe (Power-Law with Exponential Cutoff): Includes an exponential cutoff at $40 M_\odot$ to account for the PISN mass gap.
   • Parameters: Power-law index $\alpha$, mass ratio index $\beta$ (where $dN/dq \propto q^\beta$), and merger rate $R_{PL}$.
2. Hierarchical Mergers (HM):
   • Based on N-body simulations of dense stellar clusters (Ref. [57]).
   • The primary mass distribution is modeled with a scaling parameter $\lambda$ to allow flexibility (stretching/squeezing) of the mass distribution to fit data.
   • Parameters: $\lambda$ and merger rate $R_{HM}$.
3. Primordial Black Holes (PBHs):
   • Includes both monochromatic (single mass) and lognormal mass distributions.
   • Merger rates account for three channels: unperturbed early binaries, binary-single interactions inside dark matter halos, and late-time gravitational-wave captures.
   • Parameters: Characteristic mass ($m_{PBH}$ or $m_c$) and the fraction of dark matter in PBHs ($f_{PBH}$).

Statistical Analysis

• Likelihood Function: The authors use a Poisson log-likelihood to compare simulated histograms (derived from the models) against the observed LVK data bins.
• Parameter Estimation: They employ the emcee (MCMC) sampler to explore the parameter space and find the best-fit values.
• Model Comparison: They calculate the Bayesian Evidence ($Z$) and the Bayes Factor ($BF$) to compare the statistical preference of complex models (e.g., 1G + HM + PBH) against simpler baselines (e.g., 1G only). They also report the likelihood ratio ($-2\Delta \ln L$).

3. Key Contributions

1. Comprehensive Population Synthesis: The study is one of the first to simultaneously fit a three-component model (1G + HM + PBH) against the full O1–O4a LVK dataset, rather than testing populations in isolation.
2. Quantification of Hierarchical Mergers: The authors provide strong statistical evidence that a significant fraction of observed BBHs are hierarchical mergers, challenging the hypothesis that all detections are first-generation stellar remnants.
3. PBH Constraints: They derive updated constraints on the abundance of stellar-mass Primordial Black Holes ($f_{PBH}$) based on the specific mass and redshift distributions observed in GWTC-4.0.
4. Mass Gap Analysis: The paper addresses the "mass gap" issue, showing that while 1G models struggle to explain events with $m_1 > 40 M_\odot$, the inclusion of HM and PBH components naturally explains the observed peak in the $30-45 M_\odot$ range.

4. Key Results

Model Preferences

• 1G vs. 1G+HM: The data strongly prefers models including Hierarchical Mergers.
  • Comparing PLe & HM vs. PLe alone: $\ln BF \approx 71$ (decisive evidence for HM).
  • Comparing PL & HM vs. PL alone: $\ln BF \approx 66$.
  • The best-fit models suggest that 10% to 50% of the total merger rate at low redshifts is of hierarchical origin.
• Inclusion of PBHs: Adding a PBH component further improves the fit significantly.
  • PLe & HM & Lognorm PBH vs. PLe & HM: $\ln BF \approx 97$ (extremely strong evidence).
  • The likelihood ratio improves by $\Delta (-2\ln L) \approx 66$ when adding the PBH component.

Parameter Constraints

• Hierarchical Mergers: The best-fit models require a stretching of the HM mass distribution ($\lambda \approx 1.2 - 1.5$), suggesting the presence of third-generation (or higher) black holes in the sample.
• PBH Abundance:
  • The data favors a PBH population with a characteristic mass in the range $30 - 45 M_\odot$.
  • The inferred fraction of dark matter in PBHs is $f_{PBH} \simeq 3 \times 10^{-3} - 2 \times 10^{-2}$ (0.3% to 2%).
  • This corresponds to a few percent of the observed LVH events being PBH mergers.
• Mass Distribution: The lognormal PBH distribution is statistically preferred over the monochromatic distribution ($\ln BF \approx 17$), as it better matches the spread of masses observed in the $30-40 M_\odot$ peak.

Merger Rates

• The inclusion of HM and PBH populations allows the 1G merger rate ($R_{PL}$) to be lower (down to $\sim 10-20$ Gpc$^{-3}$yr$^{-1}$) compared to models where 1G must account for all events.
• The HM merger rate is estimated at $R_{HM} \approx 5 - 20$ Gpc$^{-3}$yr$^{-1}$.

5. Significance

• Origin of Massive Black Holes: The results strongly suggest that the population of black holes with masses $> 40 M_\odot$ (which are difficult to form via standard stellar collapse due to the PISN gap) are likely products of hierarchical mergers or primordial origins.
• Dark Matter Constraints: The study places competitive constraints on the abundance of stellar-mass PBHs, ruling out the possibility that they constitute 100% of dark matter in the $5-80 M_\odot$range, but confirming they could constitute a small but non-negligible fraction ($\sim 1%$).
• Future Observations: The authors emphasize that future runs (O5 and beyond) and the inclusion of spin measurements will be crucial. Since PBHs are expected to have low spins compared to hierarchical mergers (which can retain high spins), spin data will help break the degeneracy between the HM and PBH populations.
• Methodological Robustness: The findings remain robust even when varying prior assumptions (e.g., allowing negative $\beta$ values or wider rate priors), indicating that the statistical preference for a multi-population model is a genuine feature of the data rather than an artifact of modeling choices.

In conclusion, the paper argues that a single population of first-generation stellar black holes is insufficient to explain the LVK data. A combination of first-generation stars, hierarchical mergers in dense clusters, and a sub-dominant population of primordial black holes provides the most statistically consistent explanation for the observed mass and redshift distributions.