The AGORA High-resolution Galaxy Simulations Comparison Project. X: Formation and Evolution of Galaxies at the High-redshift Frontier

TL;DR

This paper compares high-resolution galaxy formation simulations at high redshift, evaluating their consistency and sensitivity to feedback and dust physics, to better understand early galaxy evolution.

Contribution

It provides a systematic comparison of multiple simulation codes for high-redshift galaxies, highlighting their convergence and differences without additional physics.

Findings

Simulations show consistent stellar mass results but differ in metallicity.

Massive halos reproduce observed properties at z=10-12 without extra physics.

Dust-to-metal ratio influences UV luminosity; dust absence increases UV brightness.

Abstract

Recent observations from JWST have revealed unexpectedly luminous galaxies, exhibiting stellar masses and luminosities significantly higher than predicted by theoretical models at Cosmic Dawn. In this study, we present a suite of cosmological zoom-in simulations targeting high-redshift ($z \geq 10$) galaxies with dark matter halo masses in the range $10^{10} - 10^{11}\ {\rm M}_{\odot}$ at $z=10$, using state-of-the-art galaxy formation simulation codes (Enzo, Ramses, Changa, Gadget-3, Gadget-4, and Gizmo). This study aims to evaluate the convergence of the participating codes and their reproducibility of high-redshift galaxies with the galaxy formation model calibrated at relatively low redshift, without additional physics for high-redshift environments. The subgrid physics follows the AGORA CosmoRun framework, with adjustments to resolution and initial conditions to emulate similar physical environments in the early universe. The participating codes show consistent results for key galaxy properties (e.g., stellar mass), but also reveal notable differences (e.g., metallicity), indicating that galaxy properties at high redshifts are highly sensitive to the feedback implementation of the simulation. Massive halos (${\rm M}_{\rm halo}\geq5\times10^{10}\,{\rm M}_{\odot}$ at $z=10$) succeed in reproducing observed stellar masses, metallicities, and UV luminosities at $10\leq z\leq12$ without requiring additional subgrid physics, but tend to underpredict those properties at higher redshift. We also find that varying the dust-to-metal ratio modestly affects UV luminosity of simulated galaxies, whereas the absence of dust significantly enhances it. In future work, higher-resolution simulations will be conducted to better understand the formation and evolution of galaxies at Cosmic Dawn.