🔭 astrophysics

Extreme Values of Black Hole to Stellar Mass Ratio for High-Redshift Galaxies

Using extreme-value statistics on JWST data for galaxies at redshifts $3.5 \lesssim z \lesssim 8.5$ , this study predicts a black hole to stellar mass ratio of approximately 0.24, a finding consistent with the highest observed values in high-redshift galaxies.

Original authors: Cameron Heather, Teeraparb Chantavat, Siri Chongchitnan, Joseph Silk

Published 2026-04-06

📖 4 min read☕ Coffee break read

CC BY 4.0

Original authors: Cameron Heather, Teeraparb Chantavat, Siri Chongchitnan, Joseph Silk

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 or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine the early universe as a bustling, chaotic construction site just a billion years after the Big Bang. For decades, astronomers believed that in this early era, the "bosses" of the galaxies (supermassive black holes) were tiny compared to the "workers" (the stars). They thought the ratio was like a giant skyscraper with a tiny ant on top—maybe 1 black hole for every 1,000 stars.

But then, the James Webb Space Telescope (JWST) arrived with its powerful new eyes. It looked back in time and found something shocking: in some of the earliest galaxies, the black holes weren't just ants; they were elephants. In some cases, the black hole was nearly as massive as all the stars in the galaxy combined.

This discovery broke the old rulebook. Scientists asked: How did these black holes grow so big, so fast, without eating all the stars?

The New Idea: It's Not Magic, It's Math

This paper by Cameron Heather and colleagues suggests we don't need to invent new laws of physics or "exotic" super-accelerators to explain this. Instead, they used a branch of math called Extreme-Value Statistics.

Here is the simple analogy:

The "Tallest Person" Analogy

Imagine you want to know how tall the tallest person in the world is.

If you look at a small town of 100 people, the tallest might be 6 feet.
If you look at a city of 1 million people, the tallest might be 7 feet.
If you look at the entire human population of 8 billion, you might find someone 8 feet tall.

You don't need to assume there is a "giant gene" that makes people grow to 8 feet. You just need to realize that when you sample a huge number of people, the extremes (the tallest, the heaviest) naturally get more extreme.

The "Galaxy Lottery"

The authors treated the early universe like a massive lottery.

The Tickets: Every galaxy is a ticket. Most tickets win a "normal" amount of stars and a "normal" black hole.
The Jackpot: But because there are so many galaxies in the early universe, the laws of probability say that somewhere, someone must have won the jackpot: a galaxy with an unusually huge number of stars AND an unusually huge black hole.
The Ratio: When you take that specific "jackpot" galaxy (the one with the biggest black hole) and divide it by the "jackpot" galaxy for stars (the one with the most stars), you get a surprisingly high ratio.

What They Did

Instead of guessing, the team used a mathematical tool (called the Weibull distribution, which sounds fancy but is just a way to predict the "worst-case" or "best-case" scenarios) to calculate what the absolute maximums should look like.

They looked at data for galaxies between redshifts 3.5 and 8.5 (very, very far back in time).
They calculated the "ceiling" for how big a black hole could be and how big a galaxy could be.
They found that the ratio of Black Hole Mass to Star Mass in these extreme cases should be around 0.24 (or roughly 1 black hole for every 4 stars).

The Result: The Math Matches the Telescope

When they compared their mathematical prediction (0.24) with the actual data JWST found, it was a perfect match.

The "elephants" (the huge black holes) that JWST saw weren't breaking the rules of physics. They were just the statistical outliers—the "tallest people" in the universe's crowd. The fact that we see them is exactly what you'd expect if you look at enough galaxies.

Why This Matters

No New Physics Needed: We don't need to invent "super-accretion" or weird new stars to explain these black holes. Standard growth, combined with the sheer number of galaxies, explains it.
The "Extreme" Connection: The paper suggests that the galaxies with the biggest black holes are also the ones with the most stars. It's a "rich get richer" scenario where the most extreme environments produce the most extreme objects.
Future Work: The authors note that if we find a black hole that is equal in mass to its stars (a ratio of 1.0) at even higher redshifts, then we might need to rethink things. But for now, the math holds up.

In a nutshell: The universe is so big and has so many galaxies that it's statistically inevitable that some of them will have black holes that look disproportionately huge. We just finally have the telescope to see them.

1. Problem Statement

Recent observations by the James Webb Space Telescope (JWST) have revealed supermassive black holes (SMBHs) in high-redshift galaxies ( $z \gtrsim 5$ ) that appear disproportionately massive compared to their host galaxies.

The Anomaly: Observed black hole-to-stellar mass ratios ( $M_{BH}/M_*$ ) reach $\sim 10\%$ at $z \gtrsim 5$ , and potentially $\sim 1$ at $z > 10$ .
The Conflict: These values starkly contrast with the local universe ratio of $\sim 10^{-3}$ (Häring & Rix 2004) and challenge standard co-evolution models.
Current Explanations: Previous attempts to explain this require "exotic" physics, such as super-Eddington accretion, top-heavy initial mass functions, or seeding by extreme Population III stars.
The Research Question: Can these extreme ratios be explained using standard statistical methods applied to empirical mass functions, without invoking new physics?

2. Methodology

The authors employ Extreme-Value Statistics (EVS), specifically the Generalized Extreme Value (GEV) formalism (also known as the block maxima method), to model the distribution of the most massive objects in a given volume.

Input Data:
- Black Hole Mass Function (BHMF): Based on A. J. Taylor et al. (2025), fitted with a Schechter function. The authors assume this function does not evolve significantly between $z \sim 4$ and $z \sim 8$ , effectively extrapolating the tail of the distribution.
- Stellar Mass Function (SMF): Based on R. Navarro-Carrera et al. (2024), also a Schechter function, corrected for Eddington bias. This function is analyzed in redshift bins from $z=3.5$ to $8.5$.
- Cosmology: Flat $\Lambda$ CDM ( $H_0 = 70$ km s $^{-1}$ Mpc $^{-1}$ , $\Omega_m = 0.32$ , $\Omega_\Lambda = 0.68$ ).
Statistical Framework:
1. Volume Division: The observational volume is divided into $N$ distinct blocks.
2. Probability Calculation: The probability $P_0$ that no galaxy in a volume $V$ exceeds a mass $M$ is modeled as a Poisson distribution: $P_0(M) = \exp(-n(>M)V)$ .
3. GEV Approximation: In the large- $N$ limit, this distribution converges to a GEV distribution:
  $G(M) = \exp(-(1+\gamma y)^{-1/\gamma})$
  where $y = (M-\alpha)/\beta$ .
4. Distribution Identification: The authors determined the shape parameter $\gamma$ for both mass functions. They found $\gamma \in [-0.15, 0]$ , indicating that the extreme distributions are well-approximated by the Weibull distribution (a specific case of GEV where $\gamma < 0$ ).
Simulation:
- The authors sampled $10^4$ instances from the derived Weibull distributions for both $M_{BH}$ and $M_*$ .
- They calculated the ratio $M_{BH}/M_*$ for each sample to generate a distribution of extreme ratios across redshift bins ( $3.5 \leq z \leq 8.5$ ).

3. Key Contributions

Statistical Explanation of High Ratios: The paper demonstrates that the observed high $M_{BH}/M_*$ ratios are a natural statistical consequence of sampling the tails of empirical mass functions, rather than necessarily requiring exotic growth mechanisms.
Quantitative Predictions: The study provides specific predicted median ratios and confidence intervals (95th and 99th percentiles) for high-redshift galaxies.
Validation against JWST: The model successfully reproduces the highest observed values from JWST (e.g., Kokorev et al. 2023; Matthee et al. 2025; Yue et al. 2024) within the 95% confidence interval.

4. Results

Predicted Ratios:
- The average predicted $M_{BH}/M_*$ ratio across the redshift range $z \sim 4-8$ is $\sim 0.24$ .
- The median ratio varies by redshift bin, ranging from 0.18 to 0.35.
- The highest predicted median occurs in the $z \in [5.5, 6.5]$ bin ( $\sim 0.35$ ), driven by non-monotonic evolution in the stellar mass function characteristic mass.
Comparison with Observations:
- The observed data points from JWST fall within the predicted 95% and 99% confidence intervals.
- Specifically, the data from Yue et al. (2024) at $z \sim 5.5-6.5$ shows a slight tension with the median prediction but remains consistent with the 95th percentile.
- The analysis of the $z \sim 11.5$ object (Bogdán et al. 2024) suggests that if confirmed, it would require a near-extreme stellar mass, potentially exceeding the standard EVS prediction derived from current Schechter fits.
Distribution Shape: The extreme mass distributions for both black holes and stars are confirmed to follow Weibull distributions, allowing for robust sampling of the "tails" of the mass functions.

5. Significance and Implications

Paradigm Shift: This work suggests that the "over-massive" black holes seen by JWST may not require new physics (like super-Eddington accretion) to explain their existence. Instead, they represent the statistical extremes of the existing population distribution.
Methodological Utility: The application of EVS provides a powerful tool to distinguish between statistical outliers and genuine physical anomalies. It establishes a baseline expectation for the most extreme objects in the early universe.
Future Directions:
- The authors note that objects at $z > 10$ with ratios $\sim 1$ might still require additional factors (e.g., heavy seeds, AGN feedback suppressing star formation) if they exceed the EVS predictions.
- The formalism can be applied to lower redshifts ( $z \sim 1$ ) to study the decoupling of black hole growth and stellar assembly, provided systematic uncertainties in mass measurements are resolved.
- Future work will explore ratios for non-extreme populations and refine mass function modeling (e.g., establishing lower cut-offs for convergence).

In conclusion, Heather et al. provide a compelling statistical framework that reconciles JWST's discovery of massive high-redshift black holes with standard cosmological models, attributing the phenomenon to the natural variance in the tails of mass functions rather than a breakdown of current physical theories.