🔭 astrophysics

Empirical estimates of how massive galaxies can be in {\Lambda}CDM

By applying Extreme Value Statistics to observed galaxy data while accounting for measurement uncertainties and scatter in the stellar-halo mass relation, this study demonstrates that the most massive galaxies observed across cosmic time remain consistent with $\Lambda$ CDM predictions, despite initial appearances of exceeding theoretical limits at high redshifts.

Original authors: Miguel Enríquez-Vargas, Aldo Rodríguez-Puebla, Aditya Manuwal, L. Y. Aaron Yung, Vladimir Avila-Reese, Carlo Cannarozzo

Published 2026-05-12

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Miguel Enríquez-Vargas, Aldo Rodríguez-Puebla, Aditya Manuwal, L. Y. Aaron Yung, Vladimir Avila-Reese, Carlo Cannarozzo

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 universe as a giant, cosmic construction site. For decades, the standard blueprint for this site—called the ΛCDM model—has told us exactly how big the buildings (galaxies) can get based on the amount of raw material (dark matter and gas) available in the neighborhood.

Recently, the James Webb Space Telescope (JWST) started taking photos of the construction site in its earliest days. It found some buildings that looked surprisingly huge for such a young neighborhood. This made some scientists worry: "Did we get the blueprint wrong? Is our theory of the universe broken?"

This paper is like a team of expert auditors coming in to check the math. They ask: "Are these buildings actually as big as they look, or are we just seeing them through a distorted lens?"

Here is what they found, explained simply:

1. The "Distorted Lens" (Measurement Errors)

When you look at a distant object, it's hard to tell exactly how heavy it is. The paper argues that our current way of weighing these ancient galaxies has a "glitch."

The Analogy: Imagine you are trying to guess the weight of people in a crowded room by looking at them from far away. Because the crowd is so dense and the view is blurry, you tend to guess that the lighter people are heavier than they really are. You accidentally "overestimate" the weight of the crowd.
The Paper's Claim: In astronomy, this is called Eddington bias. Because the number of massive galaxies drops off sharply (there are very few super-massive ones), any small error in measurement tends to push normal galaxies into the "super-massive" category. The paper shows that when you fix this "glitch" and correct for the blur, those "giant" galaxies aren't actually as massive as we thought. They shrink down to sizes that fit perfectly within the standard blueprint.

2. The "Sorting Machine" (Matching Galaxies to Halos)

The universe builds galaxies inside invisible bubbles of dark matter called halos. Think of a halo as a giant, invisible parking garage, and the galaxy is the car parked inside.

The Old Worry: Some scientists thought, "If the biggest car is bigger than the biggest parking garage, the blueprint is wrong."
The Paper's Fix: The authors used a statistical tool called Extreme Value Statistics. Think of this as a way to predict the absolute largest car you could possibly find in a specific number of parking garages.
The Result: They found that if you assume the biggest garage always holds the biggest car (a logical ranking), the math works out perfectly. The biggest galaxies fit inside the biggest halos without breaking the rules. The "tension" only appeared because we were mixing up random pairs (like putting a Ferrari in a tiny compact car garage) or because of the measurement errors mentioned above.

3. The "Fuel Tank" (Gas vs. Stars)

To build a galaxy, you need gas. The paper looked at how much gas is available in these giant halos versus how much has turned into stars.

The Analogy: Imagine a gas station with a huge tank of fuel. The question is: How much of that fuel has actually been put into cars?
The Finding:
- In the middle of the universe's history (2 to 6 billion years after the Big Bang): The galaxies were incredibly efficient. They turned almost all their available gas into stars. It was like a construction crew working at 100% speed, using up the entire fuel tank. This explains why we see such bright, massive galaxies during this specific era.
- In the very early universe (before 2 billion years): The galaxies were actually less efficient. They hadn't turned all their gas into stars yet. They were still in the "starburst" phase, building up.
- In the recent universe: The galaxies are less efficient again, with lots of fuel left over.

4. The "Dust Cloud" (Why they look so bright)

The paper also looked at how bright these galaxies are in ultraviolet light.

The Analogy: Imagine a lightbulb behind a dirty window. If you don't account for the dirt, you might think the bulb is incredibly powerful. But if you clean the window (account for dust), the bulb might just be normal.
The Result: When the authors accounted for dust and the fact that some gas is already locked up in stars, their model predicted that the brightest galaxies should look a certain way.
- For small survey areas (like looking through a keyhole), their model matched the JWST photos perfectly.
- For huge survey areas (looking at the whole sky), the model predicted galaxies that were slightly dimmer than what JWST saw. This suggests that while the mass of the galaxies fits the blueprint, there might still be some small mysteries about just how bright the biggest ones are in the largest surveys.

The Bottom Line

The paper concludes that the standard blueprint (ΛCDM) is likely still correct.

The "impossible" galaxies aren't breaking the laws of physics; they were just:

Over-weighed due to measurement errors (the distorted lens).
Mis-matched with their dark matter homes in our statistical models.

Once you clean up the data and use the right statistical tools, the most massive galaxies in the universe fit right where the theory says they should. The universe isn't broken; we just needed to adjust our glasses.

Technical Summary: Empirical estimates of how massive galaxies can be in ΛCDM

Problem Statement
Recent observations from the James Webb Space Telescope (JWST) have revealed a population of unexpectedly massive galaxies at high redshifts ( $6 \lesssim z \lesssim 10$ ), with inferred stellar masses exceeding $M_* \sim 10^{10} M_\odot$ . These findings have sparked debate regarding their compatibility with the standard $\Lambda$ CDM cosmological model. Previous studies suggested that such galaxies might require star formation efficiencies approaching or exceeding the available cosmic baryon budget, potentially challenging the hierarchical assembly of structures predicted by $\Lambda$ CDM. However, these comparisons often rely on observed stellar mass functions (GSMF) without fully accounting for systematic uncertainties, such as the Eddington bias arising from random errors in stellar mass estimates and the intrinsic scatter in the stellar-to-halo mass relation (SHMR). This paper addresses whether the apparent tension between JWST observations and $\Lambda$ CDM persists when these statistical and observational effects are rigorously incorporated.

Methodology
The authors employ an empirical framework based on Extreme Value Statistics (EVS) to estimate the probabilistic limits of the most massive galaxies and dark matter halos as a function of redshift and survey area.

EVS Framework: Instead of relying on specific galaxy formation models, the study uses the observed GSMF and the Halo Mass Function (HMF) to derive the probability distribution of the most extreme objects. The cumulative distribution function (CDF) for the maximum stellar mass, $M_{*,\text{max}}$ , is derived from the GSMF CDF, $F_*(M_*)$ , raised to the power of the number of independent galaxies, $N$ , in a given survey volume.
Correction for Uncertainties: The analysis distinguishes between three levels of GSMF:
- Observed ( $\phi_{*,\text{obs}}$ ): The raw data.
- True ( $\phi_{*,\text{true}}$ ): Deconvolved to remove random errors (Eddington bias) in stellar mass estimates.
- Intrinsic ( $\phi_{*,\text{int}}$ ): Further deconvolved to account for the intrinsic scatter in the SHMR.
  The authors model these errors using Gaussian distributions, with uncertainty ranges increasing with redshift (up to $\sim 0.6$ dex at $z \sim 7-8$ ).
Halo Constraints: The study calculates the theoretical maximum baryonic mass available in the most massive halos ( $M_{\text{bar,max}} = f_{\text{bar}} M_{\text{vir,max}}$ ) using seven different HMFs (e.g., Despali et al. 2016, Sheth & Tormen 1999) and Planck cosmological parameters.
Galaxy-Halo Connection: To compare galaxies and halos, the authors explore three correlation scenarios:
- Stochastic: Assuming the most massive galaxy and halo are independent random variables.
- Positively Correlated Scatter (PCS): Assuming a rank-order correspondence where the most massive halo hosts the most massive galaxy.
- Inverted Rank: The most massive halo hosts the least massive galaxy (used as a lower bound).
Baryonic Reservoirs: The study estimates the total cold gas mass ( $M_{\text{cold}}$ ) available for star formation by modeling cold streams and radiative cooling flows, distinguishing between hot and cold accretion modes.

Key Contributions and Results

Empirical Limits on Galaxy Mass: The authors provide EVS predictions for the most massive galaxies across various survey areas (from $10 \text{ arcmin}^2$ to full-sky). For a full-sky survey at $z \sim 0$ , the expected $5\sigma$ most massive galaxy has $M_* \sim 7 \times 10^{12} M_\odot$ , decreasing to $\sim 10^{10} M_\odot$ at $z \sim 16$ . These estimates are broadly consistent with current JWST observations when survey area is accounted for.
Impact of Eddington Bias: The study demonstrates that random errors in stellar mass estimates significantly inflate the observed abundance of massive galaxies. When correcting for these errors (moving from $\phi_{*,\text{obs}}$ to $\phi_{*,\text{true}}$ ), the inferred maximum stellar mass decreases by up to $\sim 0.4-0.5$ dex at $z \sim 5$ and nearly $\sim 1$ dex at $z \sim 16$ . Including the intrinsic SHMR scatter further reduces the maximum mass by an additional $\sim 0.1-0.2$ dex at $z \sim 5$ and $\sim 0.5-0.6$ dex at $z \sim 16$ .
Resolution of the "Impossible Galaxy" Tension: When comparing the intrinsic maximum stellar mass ( $M_{*,\text{max|int}}$ ) against the maximum available baryonic mass ( $M_{\text{bar,max}}$ ), the apparent tension largely disappears. Under the assumption of a perfect rank-order correspondence (PCS), no galaxies are found to exceed the theoretical baryonic limit ( $M_* < 0.16 \times M_{\text{vir}}$ ) at any redshift once measurement uncertainties are included.
Star Formation Efficiency (SFE):
- At intermediate redshifts ( $2 \lesssim z \lesssim 6$ ), the most massive galaxies approach the total cold gas reservoir ( $M_{*,\text{max|int}} \sim M_{\text{cold}}$ ), suggesting near-maximal star formation efficiencies ( $\epsilon \sim 1$ ) during this epoch, consistent with the peak of submillimeter/dusty star-forming galaxies.
- At higher redshifts ( $z \gtrsim 8$ ), the inferred efficiency drops significantly ( $\epsilon \ll 1$ ) when accounting for uncertainties, indicating that halos are not converting baryons with the extreme efficiency previously suggested.
UV Luminosity and Dust: The study models UV luminosities by accounting for dust attenuation and the fraction of gas already converted to stars. While the model reproduces observed UV brightness for small survey areas at high redshift, it predicts that the brightest galaxies in larger survey areas should be fainter than some current JWST detections, suggesting potential overestimation of luminosities or selection effects in the largest fields.
Cosmological Parameters: Variations in cosmological parameters (within Planck uncertainties) have a negligible impact on the EVS predictions for the most massive halos, indicating that the tension is not driven by cosmological parameter uncertainties.

Significance and Claims
The paper claims that the apparent conflict between JWST observations of massive high-redshift galaxies and the $\Lambda$ CDM model is largely an artifact of statistical and observational effects, specifically the Eddington bias and the intrinsic scatter in the SHMR. By applying a fully empirical EVS framework that rigorously accounts for these uncertainties and survey volumes, the authors conclude that current observations remain broadly consistent with $\Lambda$ CDM.

The study emphasizes that survey area is a fundamental parameter in defining the expected extremes of galaxy populations; small-area surveys are intrinsically biased toward lower mass extremes and cannot probe the full tail of the distribution. The authors assert that the "impossible galaxy" problem likely stems from systematic uncertainties in stellar mass estimates rather than a failure of the standard cosmological model or the need for new physics (such as Early Dark Energy). They further suggest that the most massive halos at $2 \lesssim z \lesssim 6$ may indeed host galaxies operating at near-maximal star formation efficiencies, while at higher redshifts, the efficiency is lower than previously claimed.

1. The "Distorted Lens" (Measurement Errors)

2. The "Sorting Machine" (Matching Galaxies to Halos)

3. The "Fuel Tank" (Gas vs. Stars)

4. The "Dust Cloud" (Why they look so bright)

The Bottom Line

More like this