On the Origin of the High Star-Formation Efficiency in Massive Galaxies at Cosmic Dawn

Z. L. Andalman, R. Teyssier, & A. Dekel

MNRAS • 2025


Abstract

Motivated by the early excess of bright galaxies seen by JWST, we run zoom-in cosmological simulations of a massive galaxy at Cosmic Dawn, in a halo of 1e11 Msol at z=9, using the hydro-gravitational code RAMSES at an effective resolution ~10 pc. We investigate physical mechanisms that enhance the star formation efficiencies (SFEs) at the high gas densities of the star-forming regions in this galaxy (~3000 H/cc, ~1e4 Msol/pc^2). Our fiducial star formation recipe uses a physically motivated, turbulence-based, multi-freefall model, avoiding ad hoc extrapolation from lower redshifts. By z=9, our simulated galaxy is a clumpy, thick, rotating disc with a high stellar mass ~3e9 Msol and high star formation rate ~50 Msol/yr. The high gas density makes supernova (SN) feedback less efficient, producing a high local SFE >10%. The global SFE is set by feedback-driven outflows and only weakly correlated with the local SFE. Photoionization heating makes SN feedback more efficient, but the integrated SFE always remains high. Intense accretion at Cosmic Dawn seeds turbulence that reduces local SFE, but this only weakly affects the global SFE. The star formation histories of our simulated galaxies are similar to observed massive galaxies at Cosmic Dawn, despite our limited resolution. We set the stage for future simulations which treat radiation self-consistently and use a higher effective resolution ~1 pc that captures the physics of star-forming clouds.


Data description

Star formation histories of simulated massive galaxies at Cosmic Dawn. Data available for the fiducial model and models lowPhot, highPhot, solTurb, varTurb. See README.txt for more information.



figure for Andalman+2025
Density-temperature distribution of gas in our suite of cosmological zoom-in simulations of a massive galaxy at Cosmic Dawn in RAMSES. The suite includes simulations with/without SNe feedback and across a range of photoionization feedback strengths. We overplot contours outlining the smallest phase space areas containing 75% of star formation events weighted by stellar mass (black) and 75% of SN events (red). In the bottom right panel, we show a schematic diagram of the gas evolution.