🔭 astrophysics

Reconstructing PTA measurements via early seeding of supermassive black holes

This paper investigates how early seeding mechanisms for supermassive black holes, specifically comparing direct collapse black holes and Dark Star collapse, influence the nanohertz gravitational wave background detected by pulsar timing arrays, finding that Dark Star seeds with a specific number density could dominate the signal while offering a pathway to constrain seed populations.

Original authors: Sohan Ghodla, Cosmin Ilie

Published 2026-07-28✓ Author reviewed ⓘ

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Sohan Ghodla, Cosmin Ilie

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine the universe as a giant, cosmic concert hall. For a long time, we thought the music playing there was mostly a quiet hum from the distant past. But recently, a group of astronomers using a network of ultra-precise cosmic clocks called Pulsar Timing Arrays (PTAs) has detected a new, low-frequency rumble—a "gravitational wave background." Think of this not as a single note, but as a chaotic, overlapping roar of sound coming from everywhere at once. The leading theory for what's making this noise is the collision of supermassive black holes, the giant monsters sitting at the centers of galaxies. When two galaxies crash, their central black holes eventually spiral together, creating ripples in space-time that we can now hear.

The big mystery this paper tackles is: where did these giant black holes come from in the first place? We know they exist today, but we also know that some of them were already huge when the universe was very young. It's like finding a fully grown oak tree in a garden that was just planted yesterday. How did they get so big, so fast? Scientists have proposed different "seed" theories—tiny black holes that formed early on and then grew. Some seeds were heavy and formed from collapsing gas clouds, while others were lighter seeds that grew from the collapse of strange, dark-matter-powered stars. This paper asks a simple but profound question: if these early seeds existed, would their descendants be the ones making all that noise we're hearing today?

The Cosmic Seed Hunt

The authors, Sohan Ghodla and Cosmin Ilie, decided to play a game of cosmic detective. They built a computer simulation to trace the life story of supermassive black holes, starting from their earliest "seeds" in the distant past and following them all the way to the present day. Their goal was to see which type of seed could produce enough black hole collisions to match the loud rumble detected by the PTAs.

They focused on two main theories for how these seeds were born. The first is the "Direct Collapse Black Hole" (DCBH). Imagine a massive cloud of gas that, instead of breaking apart into many small stars, collapses all at once into a single, heavy black hole seed. The second theory involves "Dark Stars." These are not normal stars made of burning gas; they are giant, fluffy spheres powered by the energy of dark matter particles smashing into each other. If these Dark Stars grew big enough and then ran out of fuel or got knocked around during galaxy mergers, they could collapse into black hole seeds.

The team ran their simulation to see how many of these seeds would need to exist, and how fast they would need to grow, to create the gravitational wave signal we see today. They treated the universe like a giant garden, planting seeds at different densities and watching how the resulting black holes would eventually pair up and merge.

The Surprising Winner

Here is where the story gets interesting. The simulation showed that the "Direct Collapse" seeds, while cool and heavy, are likely too rare to be the main source of the noise. It's like trying to fill a stadium with a roar using only a handful of people; even if they shout, the crowd isn't loud enough. The conditions required to form these seeds are so strict and rare that there just aren't enough of them to explain the PTA signal.

However, the "Dark Star" seeds told a different story. The authors found that if the universe was seeded with Dark Stars at a specific density—about one seed for every 1,000 cubic megaparsecs (a cubic megaparsec is a huge chunk of space, roughly 3.26 million light-years on a side)—their descendants would perfectly explain the gravitational wave background. In this scenario, the Dark Stars grow into massive black holes, which then drift into galaxies, pair up with other black holes, and spiral together. The noise from these collisions would be strong enough to dominate the signal we are hearing right now.

Setting the Limits

The paper doesn't just guess; it puts a fence around the possibilities. By comparing their simulation results with the actual data from the PTAs, the authors calculated an upper limit. They found that if there were too many of these seeds (more than about 0.01 to 0.1 per cubic megaparsec), the gravitational wave signal would be so loud that it would drown out everything else and contradict what we actually observe. So, while Dark Stars are a strong candidate for the "loud" black holes, we know they can't be too common.

In short, this research suggests that the cosmic roar we're hearing might be the echo of ancient Dark Stars collapsing and their giant black hole children fighting it out in the centers of galaxies. It's a thrilling possibility that links the invisible dark matter of the early universe to the loud, rhythmic waves we can detect today, offering a new way to understand how the giants of our cosmos came to be.

Technical Summary: Reconstructing PTA measurements via early seeding of supermassive black holes

Problem Statement
Recent detections of a stochastic gravitational wave background (SGWB) in the nanohertz (nHz) band by Pulsar Timing Arrays (PTAs) are primarily attributed to a cosmological population of inspiraling supermassive black hole (SMBH) binaries. However, observations of unexpectedly massive SMBHs at high redshifts ( $z \gtrsim 6$ ), such as CAPERS-LRD-z9 and UHZ1, challenge standard formation pathways where SMBHs grow from light seeds via accretion. These high-redshift objects suggest the existence of "heavy" seeds formed early in the Universe. While mechanisms like Direct Collapse Black Holes (DCBHs) and the collapse of Supermassive Dark Stars (SMDSs) have been proposed to explain these early massive objects, their specific contribution to the low-redshift PTA signal remains unquantified. This paper investigates whether early-seeded SMBHs, specifically those descending from SMDSs or DCBHs, can dominate the observed PTA signal and uses PTA data to constrain the number density of such seeds.

Methodology
The authors develop a formalism to calculate the SGWB spectrum ( $\Omega_{GW}$ ) produced by the inspiral and merger of SMBH descendants from early seeds. The framework involves the following steps:

Halo and Seed Population Modeling:
- The comoving number density of dark matter halos is modeled using the Extended Press–Schechter (EPS) formalism.
- The probability of a halo hosting an SMBH seed is determined by a threshold halo mass ( $M_{th}$ ). The authors assume a log-normal distribution for $M_{th}$ , truncated below a minimum mass ( $M_{min}$ ).
- Two seeding scenarios are compared:
  - SMDSs: Forming in halos with $M \sim 10^6 - 10^8 M_\odot$ at $z \approx 10-30$ , powered by WIMP annihilation.
  - DCBHs: Forming in halos with $M \gtrsim 2 \times 10^7 M_\odot$ at $z \approx 8-17$ via pristine gas collapse.
- The seed number density ( $n_{BH}$ ) is treated as a free parameter to be constrained, while the mean threshold halo mass ( $\mu_H$ ) is varied.
Growth and Evolution:
- Seeds are initialized at $z_{start} = 10$ with masses drawn from a log-normal distribution ( $10^4 - 10^6 M_\odot$ ).
- Subsequent growth is driven by accretion, calibrated to match the empirical black hole–halo mass relation observed at low redshift. The model assumes accretion dominates growth until $z \approx 1$ , after which mergers become relevant.
- Halo mergers are tracked using the EPS merger rate kernel. The authors calculate the probability that a halo merger leads to an SMBH merger, accounting for a delay time ( $\tau$ ) modeled as a log-normal distribution.
Gravitational Wave Calculation:
- The SGWB energy density spectrum $\Omega_{GW}(f)$ is calculated by integrating the merger rate density over redshift and binary mass, using phenomenological waveforms for the inspiral and merger phases.
- The calculation assumes GW-driven orbital evolution ( $\Omega_{GW} \propto f^{2/3}$ ) and neglects environmental hardening effects for the baseline prediction.

Key Contributions and Results

Dominance of SMDS-Seeded SMBHs: The study finds that if SMDSs form with a comoving number density of $n_{BH} \sim \mathcal{O}(10^{-3}) \text{ Mpc}^{-3}$ , their descendants can be the dominant contributor to the PTA signal. Specifically, binaries with total masses $m_1 + m_2 \gtrsim 10^9 M_\odot$ dominate the signal in the PTA band.
Sub-dominant Role of DCBHs: In contrast, DCBHs, with expected number densities of $\mathcal{O}(10^{-7} - 10^{-6}) \text{ Mpc}^{-3}$ , contribute only marginally to the SGWB due to their rarity and the restrictive conditions required for their formation.
Scaling Relations: The authors derive that the SGWB amplitude scales approximately as $\Omega_{GW} \propto \mu_H^{2.8} \times n_{BH}^2$ . This relationship allows for the derivation of upper bounds on seed abundance.
Upper Bounds on Seed Density: By comparing their predictions with existing PTA data (specifically the NANOGrav 15-year dataset), the authors establish an upper limit on the number density of high-redshift seeds. They conclude that $n_{BH} \lesssim \mathcal{O}(10^{-2} - 10^{-1}) \text{ Mpc}^{-3}$ . Densities significantly higher than this would overproduce the observed SGWB energy density.
Merger Frequency: The analysis indicates that for the low-number-density populations considered, SMBH mergers are infrequent. The number density of PTA-relevant SMBHs is only weakly affected by mergers because accretion-driven growth into the high-mass regime outpaces merger-driven depletion.

Significance and Claims
The paper claims to provide a method for using current PTA data to observationally constrain the number density of early-seeded SMBHs. It demonstrates that:

In the WIMP dark matter paradigm, SMDSs could plausibly seed a significant fraction of the SMBHs responsible for the PTA signal, provided their number density is within the derived bounds.
The scarcity of DCBHs likely prevents them from being the primary source of the observed nHz background.
The observed PTA signal places a strict upper bound on the abundance of heavy seeds at high redshifts ( $z \gtrsim 10$ ), offering a new observational test for early Universe seeding mechanisms.

The authors note that while their fiducial models for SMDS and DCBH abundances fall well below the derived upper limits, the framework allows for future refinement as PTA data improves and as theoretical understanding of seed formation evolves. The work does not claim to definitively prove the SMDS hypothesis but rather establishes the parameter space in which such seeds could dominate the gravitational wave background without violating current observational constraints.

Technical Summary: Reconstructing PTA measurements via early seeding of supermassive black holes

More like this