Early JWST observations revealed an excess of ultraviolet-bright, massive, and morphologically mature galaxies at redshifts z≳7z\gtrsim7–14 relative to canonical galaxy-formation predictions in the standard Λ\LambdaCDM cosmology. Shen et al. (2025) showed, via the first large-volume cosmological hydrodynamic simulations of an early dark energy (EDE) cosmology, that EDE reproduces the observed UV luminosity and stellar mass functions without recalibrating the galaxy-formation model. However, alternative baryonic prescriptions — top-heavy initial mass functions, boosted star-formation efficiency, bursty star formation, modified dust and nebular assumptions — can also enhance the predicted abundances, while the allowed EDE parameter space itself remains contested by CMB and large-scale-structure data, and both EDE and Λ\LambdaCDM runs underpredict massive quenched galaxies at z≳3z\gtrsim3 by an order of magnitude, exposing supermassive-black-hole (SMBH) feedback model limitations. The open problem is to determine, with quantified statistical confidence, whether the high-redshift galaxy population requires early-universe physics that enhances small-scale matter power or is attainable within Λ\LambdaCDM with baryonic physics, and to identify and exploit observables — galaxy clustering and bias, the reionization history and ionization topology, and Lyman-alpha forest and 21-cm small-scale power — that break the degeneracy between these two families of explanations.

Background

The concordance cosmological model (Λ\LambdaCDM: a cosmological constant Λ\Lambda governing late-time acceleration plus cold dark matter) fixes the expansion history and the primordial density-fluctuation power spectrum, from which galaxies form by hierarchical growth of dark-matter halos. The James Webb Space Telescope (JWST) found more ultraviolet (UV)-bright and massive galaxies at z≳7z\gtrsim7 than pre-JWST models predicted, and, surprisingly early, rotating gaseous and stellar disks and quenched (star-formation-suppressed) massive galaxies. A standard summary statistic is the UV luminosity function (UVLF), the comoving number density of galaxies versus absolute UV magnitude MUVM_{\rm UV}; the stellar mass function (SMF) plays the same role for stellar mass. In Λ\LambdaCDM the bright-end abundance is limited by the rarity of massive halos and by the star-formation efficiency (SFE, the fraction of a halo's baryons converted to stars), so an observed excess can be blamed on either the cosmology (more early halos) or the baryonic physics (more light per halo).

Early dark energy (EDE) is a transient additional energy component, typically an axion-like scalar field, that contributes a fraction fEDEf_{\rm EDE} of the cosmic energy budget around a critical epoch zc∼103.5z_{\rm c}\sim10^{3.5} before recombination and then dilutes away; it was originally proposed to ease the Hubble tension (the mismatch between the Hubble constant inferred from the cosmic microwave background (CMB) and measured locally). A generic consequence of EDE-like pre-recombination modifications is a reduced sound horizon compensated by a tilted primordial spectrum, which enhances the small-scale matter power spectrum and accelerates early structure formation without disturbing late-time cosmology. Whether EDE itself remains allowed is unsettled: analyses combining Planck, ACT DR6 (data release 6 of the Atacama Cosmology Telescope), SPT-3G (the third-generation South Pole Telescope CMB camera), DESI (the Dark Energy Spectroscopic Instrument) and weak-lensing (the gravitational distortion of background galaxy shapes by intervening matter) data variously find fEDEf_{\rm EDE} constrained at the 0.05–0.15 level with dataset-dependent conclusions, from preference for Λ\LambdaCDM to mild preference for EDE (McDonough et al. 2023; Poulin et al. 2025; Khalife et al. 2025).

Shen et al. (2024) first argued with an empirical (halo-occupation-type) model that EDE could jointly explain the Hubble tension and the JWST abundance puzzles. Shen et al. (2025) then ran paired (100 cMpc)3(100\,{\rm cMpc})^3 cosmological hydrodynamic simulations (2×120032\times1200^3 resolution elements, moving-mesh code with the IllustrisTNG galaxy-formation model) in Λ\LambdaCDM and in an EDE realization fitted to ACT DR4-era CMB data with fEDE(zc)=0.179f_{\rm EDE}(z_{\rm c})=0.179, using identical subgrid physics (phenomenological prescriptions for star formation, feedback, and black holes operating below the simulation's resolution scale). The EDE run matches the JWST UVLF to z≃14z\simeq14 and the SMF to z≃9z\simeq9 while the Λ\LambdaCDM run underpredicts them by ∼0.5\sim0.5–1 dex; predictions converge to Λ\LambdaCDM at z≲3z\lesssim3, and early disk emergence is enhanced. Yet the authors emphasize that significant degeneracies remain with baryonic solutions (top-heavy initial mass function (IMF), strongly boosted SFE, extreme burstiness — large temporal fluctuations in UV output at fixed halo mass — and dust/nebular-emission assumptions), that both cosmologies underpredict massive quenched galaxies at z≳3z\gtrsim3 by at least an order of magnitude (pointing to SMBH accretion/feedback subgrid limitations), and that breaking these degeneracies requires larger-volume simulations varying subgrid parameters, on-the-fly radiative transfer (RT, i.e., photon transport computed during the simulation rather than post-processed) for reionization predictions, and exploitation of distinctive observables. Three such observables stand out. Galaxy clustering: in EDE, brighter underlying matter clustering combines with reduced halo bias (galaxies of fixed luminosity occupy less rare, more abundant halos), a signature absent from pure SFE-boost scenarios. Reionization: the epoch when UV photons from the first galaxies ionized the intergalactic medium, whose history and spatial topology respond to the enhanced early galaxy abundance. Small-scale power probes: the Lyman-alpha forest (absorption by neutral hydrogen along quasar sightlines) and the 21-cm signal (hyperfine-line emission/absorption from neutral hydrogen) have mainly been used to constrain models that suppress small-scale power (e.g., warm dark matter); the opposite, enhanced-power regime relevant to EDE remains largely unexploited. The unresolved question is which family of explanations the data actually favor, and which measurements decide it.

Problem Statement

Determine whether the statistical properties of the galaxy population at z≳6z\gtrsim6–14 revealed by JWST — the bright end of the UV luminosity function, the stellar mass function and its most massive outliers, and the number density of early stellar and gaseous disks — are consistent with standard Λ\LambdaCDM combined with plausible baryonic galaxy-formation physics (star-formation efficiency, temporal burstiness, initial mass function, dust attenuation and nebular emission; "plausible" meaning that the adopted prescriptions must be compatible with the broader astrophysical calibrations and observational constraints that exist independently of the high-redshift abundance data being explained), or require pre-recombination early-universe physics that enhances small-scale matter power, with early dark energy as the representative case; and identify which observables decide the question. The abundance of massive quenched galaxies at z≳3z\gtrsim3 is not itself a discriminating statistic between the two families: both the Λ\LambdaCDM and the EDE runs underpredict it by at least an order of magnitude, a shortfall the source attributes to SMBH subgrid limitations common to both cosmologies, making it a linked open issue in its own right rather than evidence for either family. The source's future-work discussion identifies three linked directions that would advance a solution (candidate routes, not mandatory deliverables):

(1) Simulation-based prediction with controlled degeneracies. Cosmological hydrodynamic simulation suites (of sufficient volume and resolution; possible extensions include on-the-fly radiative transfer for reionization-era predictions and zoom-in runs resolving the multiphase interstellar medium), run in Λ\LambdaCDM and in enhanced-small-scale-power cosmologies — with EDE parameters propagated from current CMB and large-scale-structure constraints, e.g., fEDEf_{\rm EDE} and nsn_{\rm s} — while systematically varying the galaxy-formation subgrid prescriptions, so that the covariance between cosmological and baryonic variations is quantified for the predicted observables entering the comparison.

(2) Discriminating observables. Quantitative predictions for statistics that differ between the families even when both match the abundance data: galaxy clustering and bias, where EDE predicts enhanced underlying matter clustering combined with reduced halo bias at fixed galaxy properties; the reionization history and ionization topology; and the imprint of enhanced small-scale power on the Lyman-alpha forest and the 21-cm signal — a regime these probes have largely not been used to constrain.

(3) Statistical discrimination on data. A joint statistical comparison with current and foreseeable observations, with quantified significance and with the dominant systematics (cosmic/sample variance, stellar-mass and photometric systematics including Eddington bias (the spurious enhancement of inferred abundances when measurement scatter pushes sources across a steep selection boundary), foreground contamination) propagated, that either demonstrates that Λ\LambdaCDM plus baryonic physics suffices at the precision of the data, establishes the requirement for enhanced small-scale power, or quantitatively demonstrates that the available data cannot yet separate the two families.

Completion is a defensible, quantitative discrimination: on the basis of a statistical comparison with the dominant systematics propagated, the answer must state which explanation the data favor, at what confidence, over which observables, and why the competing family is disfavored — or, if the data cannot yet separate the two families, demonstrate that quantitatively and identify which measurements, with what projected sensitivity, would decide the question.

References

  1. Shen, X., Zier, O., Vogelsberger, M., Boylan-Kolchin, M., Hernquist, L., Tacchella, S., Naidu, R. P. (2025). The Cosmic Rush Hour: Rapid Formation of Bright, Massive, Disky, Star-Forming Galaxies as Signatures of Early-Universe Physics. arXiv:2509.19427; doi:10.1093/mnras/stag1227 (published in Mon. Not. R. Astron. Soc. 550, issue 2, article stag1227, 2026). https://arxiv.org/abs/2509.19427
  2. Shen, X., Vogelsberger, M., Boylan-Kolchin, M., Tacchella, S., Naidu, R. P. (2024). Early Galaxies and Early Dark Energy: A Unified Solution to the Hubble Tension and Puzzles of Massive Bright Galaxies revealed by JWST. arXiv:2406.15548; doi:10.1093/mnras/stae1932 (published in Mon. Not. R. Astron. Soc. 533, 2024). https://arxiv.org/abs/2406.15548
  3. Poulin, V., Smith, T. L., Calderón, R., Simon, T. (2025). Impact of ACT DR6 and DESI DR2 for Early Dark Energy and the Hubble tension. arXiv:2505.08051; doi:10.1103/bx25-1g5d (published in Phys. Rev. D 113, 2026). https://arxiv.org/abs/2505.08051
  4. McDonough, E., Hill, J. C., Ivanov, M. M., La Posta, A., Toomey, M. W. (2023). Observational constraints on early dark energy. arXiv:2310.19899; doi:10.1142/S0218271824300039 (published in Int. J. Mod. Phys. D 33, 2024). https://arxiv.org/abs/2310.19899
  5. Khalife, A. R., Balkenhol, L., Camphuis, E., et al. (SPT-3G Collaboration) (2025). SPT-3G D1: Axion Early Dark Energy with CMB experiments and DESI observations. arXiv:2507.23355; doi:10.1103/8jjr-7hpb (published in Phys. Rev. D 113, 2026). https://arxiv.org/abs/2507.23355
  6. Comini, L., Vagnozzi, S., Loeb, A. (2026). Dark energy, spatial curvature, and star formation efficiency from JWST photometric and spectroscopic high-redshift galaxies. arXiv:2604.13866; doi:10.1016/j.jheap.2026.100626 (published in Journal of High Energy Astrophysics 53, 2026). https://arxiv.org/abs/2604.13866
  7. Weibel, A., Jespersen, C. K., Oesch, P. A., Williams, C. C., Bezanson, R., Brammer, G., Cloonan, A. P., Dayal, P., Hutter, A., Ji, Z., Maseda, M. V., Shuntov, M., Whitaker, K. E. (2025). Exploring Cosmic Dawn with PANORAMIC II: Cosmic Variance and Galaxy Clustering at z~10. arXiv:2512.14212. https://arxiv.org/abs/2512.14212
  8. Dalmasso, N., Ferrami, G., Leethochawalit, N., Ventura, E. M., Trenti, M. (2026). Accelerated evolution of galaxy host halo masses during Cosmic Dawn from deep JWST clustering. arXiv:2601.01697. https://arxiv.org/abs/2601.01697
  9. Matteri, A., Pallottini, A., Ferrara, A. (2026). Clustering constraints on super-early galaxy formation scenarios. arXiv:2605.11064. https://arxiv.org/abs/2605.11064
  10. Hodgson, S., Wells, P., Heitmann, K., Krishnan, N. (2026). Investigating The Effects of Early Dark Energy on Large-scale Structure Within the EDENS Suite. arXiv:2608.17212. https://arxiv.org/abs/2608.17212
  11. Zier, O., Smith, A., Shen, X., Liu, R., Kannan, R., Koehler, S. M., Springel, V., Pakmor, R., Vogelsberger, M., Bulichi, T.-E., Hernquist, L. (2026). Introducing the Lumina project: large-volume radiation-hydrodynamic simulations of the epochs of hydrogen and helium reionization. arXiv:2605.15310. https://arxiv.org/abs/2605.15310
  12. Samuel, J., Boylan-Kolchin, M., Feldmann, R., Hopkins, P., et al. (2026). Resolving galaxy formation in the early Universe with BonFIRE and CampFIRE. arXiv:2605.24104. https://arxiv.org/abs/2605.24104
  13. Liu, W., Wang, X., Zhan, H., Glazebrook, K., Tang, M., Trenti, M. (2026). A Portrait of the Cosmic Reionisation History in the Context of the Early Dark Energy Model. arXiv:2601.11068; doi:10.3847/1538-4357/ae3d0c (published in Astrophys. J. 999, 2026). https://arxiv.org/abs/2601.11068
  14. Lagos, C. del P., Valentino, F., Wright, R. J., de Graaff, A., Glazebrook, K., et al. (2024). The diverse star formation histories of early massive, quenched galaxies in modern galaxy formation simulations. arXiv:2409.16916. https://arxiv.org/abs/2409.16916
  15. Ji, Z., Williams, C. C., Behroozi, P., Weibel, A., Jespersen, C. K., Oesch, P. A., et al. (2026). PANORAMIC: The Dawn of Massive Quiescent Galaxies I. Number Density and Cosmic Variance from 1000 arcmin^2 NIRCam Imaging. arXiv:2604.05022. https://arxiv.org/abs/2604.05022
  16. Adi, T., Flitter, J., Kovetz, E. D. (2024). Early Dark Energy Effects on the 21-cm Signal. arXiv:2410.22424; doi:10.1103/PhysRevD.111.043515 (published in Phys. Rev. D 111, 043515, 2025). https://arxiv.org/abs/2410.22424
  17. Mazengo, D. T., Popesso, P., Marini, I., et al. (2026). What's Missing in AGN Feedback? Lessons learnt from Magneticum, IllustrisTNG and Simba. arXiv:2607.21140. https://arxiv.org/abs/2607.21140
  18. Krishnan, J. R., Abazajian, K. N. (2026). Statistics Meet Systematics: Resolution of the Massive Early JWST Galaxy Tension. arXiv:2511.13708; doi:10.1103/gjf5-3r89 (published in Phys. Rev. D 113, 2026). https://arxiv.org/abs/2511.13708
  19. Enríquez-Vargas, M., Rodríguez-Puebla, A., Manuwal, A., Yung, L. Y. A., Avila-Reese, V., Cannarozzo, C. (2026). Empirical estimates of how massive galaxies can be in ΛCDM. arXiv:2605.08353. https://arxiv.org/abs/2605.08353
  20. Shen, X., Zier, O., Smith, A., Liu, R., Kannan, R., Bulichi, T.-E., Koehler, S. M., Springel, V., Vogelsberger, M., Hernquist, L., Naidu, R. P., et al. (2026). The Lumina Project: The Demographics of Active Galactic Nuclei from Quasars to Little Red Dots at z>=3. arXiv:2605.24112. https://arxiv.org/abs/2605.24112
  21. Bhattacharjee, A., Chakraborty, A., Das, S., Maharana, A., Parashari, P. (2026). Transient Early Dark Energy-Like Dynamics as a Mechanism for Enhanced Early Structure Formation in the JWST Era. arXiv:2608.20288. https://arxiv.org/abs/2608.20288